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A simplified formula for light-deflection in the quadrupole field of moving massive 
bodies has been obtained in [ll-0], which will be applied for Gaia data reduction. 
So far, in Gaia data reduction it has been assumed that the positions of the giant 
planets should be computed at the retarded instant of time. The problem of light- 
deflection due to quadrupole field of moving planets has been re-considered in Q]. 
According to their solution, the position and velocity of the massive body have to 
be taken at retarded time. We show that the solution given in coincides with our 
simplified quadrupole formula obtained in [ll-0]. This coincidence implies that the 
positions of giant planets have, in fact, to be taken at retarded time. 
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I. INTRODUCTION 



An essential experiment of Gaia mission to test relativity concerns the new effect of 
light-deflection at giant planets due to their quadrupole field. Analytical solutions of light 
deflection by a quadrupole field of a massive body are well-known and have been investi- 
gated by many authors @-[l2|. Because these formulas are rather complicated, they imply 



massive computations of quadrupole light deflection and are too time-consuming for Gaia 
data reduction. Therefore, we have derived a simplified quadrupole formula in [l|-l3|, which 
is suitable for a time-efficient computation of quadrupole light-defiection on microarcsecond 
level of accuracy. 

In reality, while the light-signal is being emitted at a position Xq at time moment to and 
received at position Xi at a time moment ti, the massive body chances the position from 
XA{to) to XAiti). Therefore, it is not obvious at which coordinate time to < t < ti the 
position of the massive body XA{t) has to be chosen in the formula of quadrupole light- 
defiection. Up to now, in Gaia data reduction it is tacitly and implicitly assumed that the 
positions of multipoles (Jupiter and Saturn) should be computed at the retarded moments 
of time, given by the implicit light-cone equation: 



^ ^^-^ . (1) 

where Xi is the position of observer at observation time ti and x^. (s^) is the coordinate of 
massive body at retarded time moment s^. However, no theoretical proof of this assumption 
has been given. This problem has recently been solved by 0]. In our report we show that the 
solutions in and in agree with each other, implying that in the quadrupole formula 
the position of massive body has indeed to be taken at the retarded instant of time. 



II. SIMPLIFIED QUADRUPOLE FORMULA 

In a simplified quadrupole formula has been derived which takes into account only 
those terms relevant for microarcsecond level of accuracy. In this Section we will give the 
main steps and results. Consider a gravitational field of massive bodies A, and the 
positions of these individual massive bodies are xa- A light-ray is being emitted at a 
position Xq at time moment to and received at position Xi at a time moment ti. The unit 

x(ti) 

coordinate direction of the light propagation at the moment of observation reads n = - — - — -- 

\x{ti)\ 

x(t) 

and the unit tangent vector of light path at infinitely past is cr = lim . Then, in post- 

Newtonian order, the transformation cr to n is given by [9|, llO 



n 



J2s<T. + (c^^) , (2) 



where the sum runs over individual terms of various physical origin, that means monopole 
gravitational field, quadrupole field and higher multipole fields. The spherical symmetric 
part (monopole field) due to one massive body A is given, for instance, by Eq. (102) in [isj . 

Here, we are only interested at the quadrupole light deflection term. For one massive 
body A, the quadrupole light-deflection is given as follows (see Eq. (40) in ^ or Eq. (9) in 
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131): 

)2 



5o-Q = £^ (2 + 3 cos ^ - cos=^ ^) 

C 



X [(l - (cr ■ 63)^ - 4 (n^ • 63)^) + 2 (n^ ■ 63) 63 - 2 (cr ■ 63) (ua ■ 63) cr] . (3) 

Here, Ma is the mass of body A, c is the speed of hght, G is the gravitational constant, and 
the impact vector 

cLa = (T X {rf X a) , (4) 

where = a?i — xa, and the absolute value (Ia = \dA\- Furthermore, the angle -0 = 
6 (cr,r^), the unit vector along the axis of symmetry (rotational axis of massive body) is 
denoted by 63, Pa denotes the equatorial radius, and is the coefficient of second zonal 
harmonic of the gravitational field of massive body A. The unit vector is defined by 

Since the massive bodies move, the question arises at which instant of time the positions of 
the massive bodies have to be chosen. This problem has been solved in 0], and the main 
results of this work will be subject of the next Section. 



III. QUADRUPOLE FORMULA BY KOPEIKIN & MAKAROV 

A. Description of the approach 

The light-defiection at moving monopoles and quadrupoles has been re-investigated by in 
In this Subsection we describe the basic steps of this approach. The gravitational field 
is described by gap = i]ai3 + where hap is the metric perturbation in post- Minkowski 
apprximation. Using harmonic gauge, the linearized Einstein equations for the pertubation 
ha/3 are homogeneus wave equations, 

+ = 0. (6) 



A general solution of (|6]) is given in terms of multipole expansion [IJ, ll5|], where the terms in 
the perturbation hap are taken at the retarded instant of time sf. The retardation follows 
directly from the retarded (causal) solution of the homogeneus wave equations (Q. 

In the approach 0] the monopole, dipole and quadrupole terms of the general multipole 
expansion of hap, given in Eqs. (12) - (14) in are taken into account. The center of 
coordinate system is shifted from the mass center of the massive body by a spatial distance. 
Then, by means of parallel axis theorem (Huygens-Steiner theorem), they apply general 
expressions for the quadrupole moment and their time derivative, see Eq. (17) and Eq. (20) 
in g. 

Then, in yi the geodesic equation is rewritten into a considerably simplier form, given by 
Eq. (19) in [l6|. The integration leads to the expressions (28) and (29) in 14]. Using their 
developed integration method, presented in part by Eqs. (30) and (31) in they succeed 
to integrate analytically the differential equation (29) in Furthermore, in j3] all those 
terms are neglected which they proof to contribute less jhan 1 microarcsecond. The result 
of their approach is finally given by Eqs. (39) - (41) in 
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B. The analytical solution for quadrupole light-deflection 

In Eq. (44) in , the following form for the light-deflecetion of sources at infinite distances 
due to one moving quadrupole A has been given: 



X [[(eg ■ ua) - (63 ■ nriA) ] - 2 (eg ■ ua) (es ■ m^) niA] ■ (7) 

It is essential to noticed, that all time-dependent quantities in Eq. ([7]), that is xa, va and 
63, and therefore also ua and tua, are computed at retarded instant of time given by 
Eq. ([1]). It should be noticed that in Eq. ([7]) we have omitted a term proportional to the 
displacement of the planetary center from the origin of the coordinate system which is an 
artificial term of the approach and not relevant for Gaia data reduction. Furthermore, the 
unit vector tha is defined by 

iriA = cr X nA- (8) 

In Eq. ([7]) we have used our notational conventions in [2I, 0] . As mentioned in the introduc- 
tionary Section, we have replaced the vector k by cr, since cr = k + O (m), that means such a 
replacement would cause effects of higher order beyond the post-Newtonian approximation. 



C. Comparison of the formula of quadrupole light-deflection 

The quadrupole light-defiection ([7]) can be estimated by 

IKI.^(i-=^)^f . (.) 

Accordingly, for Jupiter and Saturn we have 

Jupiter 



< (1 _ ^-^J"piter -) 240 /iaS, (10) 



Saturn \ ^ f ^ ^ ' ''^Saturn \ 



< {1- ™° j95/xas, (11) 



while for other massive bodies of the solar system we find considerably smaller values. 
Furthermore, the orbital speed of these planets in respect to the barycenter is of the order 
Va ~ 10~^ c. Hence, the contribution of the velocity-term in the quadrupole light-defiection 
is by far less than about 0.1 ^as and can be neglected for Gaia astrometric accuracy. Thus, 
the quadrupole light-defiection formula ([7]) can be simplified as follows: 

x^A A GMa a Pa 

X [[(63 ■ Ua^ - (63 ■ rriAf] nA-2 {e^ ■ ua) (eg ■ m^) uia] ■ (12) 

Eq. f lT^ determines the quadrupole light-defiection for moving massive bodies, while the 
rotational axis 63 and the unit vectors ua and tua have to be computed at retarded instant 
of time given by Eq. ([T]) . 
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The expression given by Eq. ( IT^ in coincides with our expression given in Eq. , see 
This can be shown by means of the relation cos^ x = 1 — sin^ x, that is (63 ■ m^)^ = 
1 — (63 X m^i)^. Then, by taking into account 63 x = 63 x (cr x n^) we obtain the 
relation 

(63 ■ ruAf = 1 - (63 ■ uaY - (63 ■ o-f . (13) 



Inserting ( IT^ into (IT^ we obtain 

Ma jaQ 

C2 ^3 



OCTo = -4 — Jn — r 



X [[1 - 2 (63 ■ tia) - (eg ■ cr) ] riA + 2 (eg ■ tt-a) (eg ■ /tia) ttia] • (14) 

Furthermore, for the vector e3 we have to use the linear combination in terms of the unit 
vectors d^, cr and tua, that means eg = (63 ■ cr) cr + (63 ■ Ua) tia + (eg ■ tua) "TUa, from 
which we conclude 

(eg ■ niA) rriA = eg - (eg ■ cr) cr - (eg • Ua) Ua ■ (15) 
Inserting f lTSj) into f|T^ . we obtain 

r^A _ 4GMA A^l 

X [(1 - (cr ■ eg)^ - 4 (nA ■ eg)^) riA + 2 (ua ■ eg) eg - 2 (cr • eg) {ua ■ eg) cr] . (16) 

This expression coincides with our quadrupole light-deflection for a massive body, given by 
Eq. 03]) if we approximate 2 + 3 cos ip — cos^ ip hy a, factor 4. This approximation follows 

d ( 
from sin?/^ = — , that means 2 + 3 cos'i/' — 0,0^ ip = 4 + 01 — ^ | . Thus, we have shown 



the coincidence of our quadrupole formula given by Eq. ([3]) (that means Eq. (40) in [2| or 
Eq. (9) in (3]) with the quadrupole formula in j^] given by Eq. f|T2|) . This coincidence implies 
that for Gaia data reduction the position of giant planets has, in fact, to be taken at the 
retarded instant of time, given by the implicit light-cone equation ([T]). 



IV. SUMMARY 

A simplified quadrupole formula for massive bodies has been derived in (1M3I , which will be 
used for a time-efficient computation of quadrupole light-deflection in Gaia data reduction. 
However, since the massive bodies move during the light-signal travelling, it is not obvious 
at which coordinate time the position of the massive body XAif) has to be taken in this 
formula of quadrupole light-defiection. So far, in Gaia data reduction it has been tacitly 
assumed that the positions of the giant planets should be computed at the retarded instant 
of time Xa{,Sy^. However, no theoretical proof of this assumption has been performed. 

Therefore, in this report, we have re-considered the results of [3] where the light-deflection 
at moving quadrupoles has been determined. We have shown that the quadrupole formula 
obtained in [3] given by Eq. f[T^ coincides with our simplified quadrupole formula obtained 
in [l|-[^ and given by Eq. ([3]). This coincidence implies that the position of massive body 
has, in fact, to be taken at the retarded instant of time. 
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